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Abstract 

, We show that the current accelerated expansion of the Universe can be explained without resort- 

ed ' 

ly-^ ' ing to dark energy. Models of generalized modified gravity, with inverse powers of the curvature 

■ can have late time accelerating attractors without conflicting with solar system experiments. We 

^ \ 

. have solved the Friedman equations for the full dynamical range of the evolution of the Universe. 

O 

fH \ This allows us to perform a detailed analysis of Supernovae data in the context of such models 

I ' that results in an excellent fit. Hence, inverse curvature gravity models represent an example of 

H 



phenomenologically viable models in which the current acceleration of the Universe is driven by 
curvature instead of dark energy. If we further include constraints on the current expansion rate 
of the Universe from the Hubble Space Telescope and on the age of the Universe from globular 



C/3 
> 

■ clusters, we obtain that the matter content of the Universe is 0.07 < < 0.21 (95% Confidence). 

Hence the inverse curvature gravity models considered can not explain the dynamics of the Universe 
just with a baryonic matter component. 

PACS:04.50.+h, 95.36.+X 
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It is now widely accepted that recent Supernovae (SNe) observations imply that our 
Universe is currently experiencing a phase of accelerated expansion This seems to be 
independently;^ confirmed by observations of clusters of galaxies 0] and the cosmic microwave 
background 3]. The accelerated expansion is usually explained through violations of the 
strong energy condition by introducing an extra component in the Einstein equations in the 
form of dark energy with an equation of state w < —1/3. However, such an explanation is 
plagued with theoretical and phenomenological problems, such as the extreme fine tuning 
of initial conditions and the so called coincidence problem P| and it is therefore natural 
to seek alternatives to dark energy as the source of the acceleration. One possibility is an 
inhomogeneous Universe with only local acceleration, albeit it is hard to explain natural 
boundary conditions for such a local void [5]. The other, that we will elaborate on in this 



Letter, is rnodjfications of gravity that turn on only at very large distances ^6^ or small 



therefore giving a geometrical origin to the accelerated expansion of the 



curvatures 
Universe. 

It was shown in [3] that a simple modification of the gravitational action adding inverse 
of curvature invariants to the Einstein-Hilbert term would naturally have effects only at 
low curvatures and therefore at late cosmological times. The simplest of such modifications 
includes just one single inverse of the curvature scalar fi'^/R, with /i a parameter with 
dimensions of mass. This results in a model governed by the Einstein-Hilbert term, i. e. usual 
gravity, for curvatures fi^ but can lead to an accelerated expansion at curvatures R < 
/i^. This simple model is equivalent to a Brans-Dicke theory ^j. Based on this equivalence 
it was subsequently proven by a number of authors that the model is in confiict with solar 
system data jg] and is unstable when matter is introduced jl^. This conclusion naturally 
extends to generalizations of this action where the Einstein-Hilbert term is supplemented 
with an arbitrary function of R, except for particular cases that could still lead to viable 
models 

With this restriction in mind, the authors of discussed a more general modification 
of gravity based on the following gravitational action 



^ = 16^ / ^'^^ 



,,4n+2 

R- ^ 



{aR^ + hP + cQY 



(1) 



where P = R^^R^^ , Q = R^^p^R^^P" , G Newton's constant. Cm the matter Lagrangian and 



g the determinant of the metric. 

In this generahzed case the equivalence with a Brans-Dicke theory is not clear and a 
more detailed analysis of modifications of Newton's potential has to be done to compare 
with solar system data. The authors of computed the corrections to Newton's law in 
these models as a perturbation around Schwarzschild geometry and found that as long as we 
include inverse powers of the Riemann tensor (c 7^ 0), Newton's law is not modified in the 
solar system at distances shorter than Vc ~ 10 pc and therefore all solar system experiments 
are well under control. Note that, as long as the Riemann tensor is present, this result is 
independent of whether we include or not inverse powers of the scalar curvature or the Ricci 
tensor squared, as they vanish in the background solution. This important result restricts 
the parameter space of phenomenologically relevant inverse curvature gravity models to the 
ones with inverse powers of the Riemann tensor squared present. Other constraints come 
from the absence of ghosts in the spectrum, requiring specific relations between b and c Q| . 
Finally we restrict our analysis in this Letter to models with n = 1. 

Let us turn now to the cosmology of models governed by the gravitational action (P). 
Assuming a cosmological setup with a spatially flat Friedmann-Robertson- Walker metric, 
ds'^ = —dt^ + a{t)da^, all models with n = 1 can be characterized by just three parameters, 
a, fi and a, given in terms of the parameters in Eq. (P) by 

_ 12a + 46 + 4c 
" " 12a + 36 + 2c' ^ ^ 

/t =/i/|12a + 36 + 2c|^/^ (3) 
a =sign(12a + 36 + 2c). (4) 

In order to write the corresponding Friedmann equation in the simplest possible way we will 
use logarithmic variables, u = \n{H / fi) and = Ina, where as usual H = ^, with a dot 
denoting the time derivative. The generalized Friedmann equation in these variables reads 

u"Viiu') + + 18a(P3(M'))'e'"(e'^"""^ - 1) = 0, (5) 

where a prime denotes the derivative with respect to N and we have defined the following 
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polynomials, 

Vi (y) = 6a^y^ + 24ay + 32 - 8a, (6) 
V^(y) = 15aV + 2a(50 - 3a)y^ 

+ 4(40 + lla)?/2 + 24(8-a)2/ + 32, (7) 
Vsiy) = ay^ + Ay + 4 . (8) 

The source is u = In [tu^ exp(— 4A^) + LJmexp{—3N)] /2, where we have defined the appropri- 
ately normalized values of the energy densities today as 

SirG Pr^mO _ - /qv, 
o p 

with p°j^ the present densities in matter and radiation and we have exploited the fact that 
the energy-momentum tensor is still covariantly conserved. This means that the source in 
Eq. (0) corresponds to the standard one with no dark energy. 

The new Friedmann equation is no longer algebraic but a second order non-linear dif- 
ferential equation. Furthermore, it becomes non-autonomous in the presence of sources, 
making its dynamical study a formidable problem. The asymptotic behavior of the system 
in vacuum was carefully studied in 3|, where it was found that, depending on the value 
of a, but irrespective of a, the system has a number of attractors, including sometimes 
singularities. The same attractor and singular points are relevant when sources are present. 
In that case however, both the value of a and the fact that the Universe is in a matter 
dominated era before the new corrections become relevant are crucial to determine the fate 
of the Universe. 

A careful analysis of the dynamical behavior of the system reveals that physically valid 
solutions only exist for certain combinations of a and a. In order to classify the different 
regions, we define the following special values of a: ai = 8/9, 0^2 = 4(11 — vT3)/27 ~ 1.095 
and as = 20(2 - V3)/3 ^ 1.786. For a < ai both signs of a result in an acceptable 
(non-singular) dynamical evolution, but nevertheless in a bad fit to Supernovae data. For 
ai < a < a2 only a = —1 leads to an acceptable expansion history, since for a = +1 
a singular point is violently approached in the past. For a2 < a the singular point is 
approached for a = —1, hence a = +1 is the only physically valid solution. In this latter 
case, when a2 < a < a^, the system goes to a stable attractor that is decelerated, thus giving 
a bad fit to SNe data, for a < 32/21 and gets accelerated for larger a. For < a there 



is no longer a stable attractor and the system smoothly goes to a singularity in the future. 
That singularity occurs earlier as a increases so that there is a limiting function ai{ujm), 
at which the singularity is reached today. It is important to stress that this singularity is 
approached in a very smooth fashion, allowing for a phenomenologically viable behavior of 
the system, as opposed to the evolution when the wrong value of a is chosen, where the 
singularity is hit almost instantaneously. Finally, for values of a > 24.9, there are stable 
attractors again but these are never accelerated and the resulting fit to SNe data is not 
acceptable. To summarize, there are two regions that give a dynamical evolution of the 
system compatible with SNe data, the low region with q;i < a < q;2, for which a — —1, and 
the high region where a2 < Q < o^, for which a = +1. 

As we have emphasized it is extremely difficult to solve the dynamics of the system 
analytically. To overcome this limitation we have performed a comprehensive numerical 
study of the model resulting in the general behavior we have outlined above. To make things 
more complicated the new Priedmann equation is extremely stiff, due to the exponentials 
in the last term. This stiffness is directly hnked to the nature of the corrections that are 
negligibly small in the far past, where the curvature is much smaller than the scale fx^. It also 
makes it essentially impossible to numerically integrate it from a radiation dominated era all 
the way to the present. In order to circumvent this problem, we have matched a perturbative 
analytical solution that tracks the solution in standard Einstein gravity in the far past to 
the corresponding numerical one in the region > 5, where the analytical solution is still 
an extremely good approximation, and the numerical codes can cope with the integration. 
Although the matching at this point is accurate below the 1% level, we emphasize that 
it is safely above the redshift range probed by SNe. The approximate solution from the 
perturbation analysis, for o; 7^ 8/9, is given by 



This is an extremely accurate solution to the full non-linear equation as long as z > 5, 
regardless of the values of a and uJm- At the boundaries between regions with different 
dynamical behavior (including ai) the sensitivity to initial conditions is large and therefore 
nothing conclusive can be said at these points. The question of sensitivity to initial condi- 
tions is a relevant one due to the non-linear nature of Priedmann equation. However due to 
the complication of any analytical study for non-negligible sources alluded to above we will 




(10) 
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defer its study to a future publication. In the present Letter we will contempt ourselves with 
the particular solution in Eq. (fTUI) that we are guaranteed tracks the standard behavior in 
Einstein gravity in the past. We further explicitly confirmed, by a numerical analysis, that 
our conclusions are not sensitive to the exact position of the matching point in the past. 

Once we have solved for the Hubble parameter as a function of the scale factor, we perform 
a fit to SNe data to get the allowed values of the different parameters defining our model. In 
principle there is a total of five parameters defining our Universe in this framework, namely 
the three parameters defining the model, fi, a and a, and the two parameters determining 
the sources, iDr,m- The absolute value of the CMB temperature, however, fixes the total 
radiation content of the Universe, constraining a),./!^. For relevant values of jj, this constraint 
makes radiation irrelevant in the analysis of SNe data. Since the intrinsic magnitude of SNe 
is a nuisance parameter in our analysis, it is not possible to determine fi as an independent 
parameter with SNe only. For a standard ACDM Universe this corresponds to the inability 
of SNe data to independently determine the Hubble constant Hq. However, we will be able 
to determine the value of a once we impose other constraints, like the measurement of the 

_ _ n 

Hubble constant by the Hubble Key Project, ifo = 72 ± 8 Km s ^ Mpc ^ Hence, this 
leaves us with just three parameters, a, a and relevant for the analysis of SNe data and 
an additional nuisance parameter in terms of the intrinsic magnitude. 

The fits are performed using the recent gold SNe data set from the last reference in 
The apparent magnitude is given by m(z) = M. + SlogP^ where = M — Slog/i + 25 
and Vl = fidi with di oc J H^^{z)dz. Note that the parameter jj, appears in the definition 
of the magnitude compared to the usual definition involving Hq The important point 
is that 'Dl{z) can now be computed solely in terms of Um and a, where jj, and the intrinsic 
magnitude have been absorbed into the nuisance parameter Ad that can be marginalized 
analytically in the probability function. 

We have performed independent two parameter fits to SNe data for each of the low and 
high regions. This results in the 1- and 2-a joint likelihoods shown in FigQ with best fit 
values given by 

low: a = 0.9, tUm = 0.105, = 184.9, (11) 
high : a = 2.15, tOm = 0.085, = 185.2. (12) 

For comparison purposes, we have also performed the fit using the standard ACDM model 
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FIG. 1: 1 and 2-a joint likelihoods on u)m and a. In the low region a = —1 whereas in the high 
region a = +1. The shaded area on the right determines the region a > 04 that is excluded because 
of a singularity being hit in the past. The diamonds denote the maximum likelihood points. 

for a spatially flat Universe and absorbing Hq as a nuisance parameter into M., resulting 
in = 184.9 for 156 data points. We further show in Fig. ^the points a^. The shaded 
area on the right side, which is bordered by a dotted line, is the exclusion zone given by 
04(0)111). Note that the contours have a sharp cut-off at ai, and 0:4. However at 03 there 
is no singularity hit violently and the 2 — 0" contour of the high region extends below 03. In 
the low region we obtain uj^ = 0.122 ± 0.034 after marginalization over a and in the high 
region Co^ = 0.075 ± 0.031. Note that our best fit points in both regions are close to the 
borders of the allowed region. This is because within the regions there is a smooth behavior 
of the likelihood, and only the dynamics of the system cuts off the likelihood space if certain 
parameter values are reached. 

If we additionally apply the HST measurement of Hq, we can determine fi and the 
matter content = l^m^^, with Hq = 100/ikm/s/Mpc. Finally we can restrict the allowed 
region in — fi a little bit more by imposing a prior on the age of the Universe with a 
mean of to = 13.4 Gyrs and a 95% confidence lower limit of 11.2 Gyrs ^] . In Fig. El we 
show the joint 1- and 2-a likelihoods in the cUm-A plane, with both priors imposed (solid 
line) and without imposing the age of the Universe prior (dashed line). On the left for the 
low region and on the right for the high region. First we recognize that fi is roughly twice 
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FIG. 2: The 1- and 2-a joint likelihoods in the — A plane, when additional priors on Hq and the 
age of the Universe are imposed. On the left for the low region and on the right for the high region. 
Diamonds are the maximum likelihoods. Further, the dashed contours are the joint likelihoods if 
we impose only the Hq prior. 



the size of the Hubble constant Ho. If we further marginalize over fi the physical matter 
content in the Universe is uj^ = 0.14 ± 0.03 and = 0.14 ± 0.04 in the low and high 
regions, respectively. Note that the matter content in the budget of the Universe is clearly 
higher than the measured baryonic content. Overall we find 0.07 < uj^ < 0.21 at the 95% 
confidence level. If we compare this number with the results from Big Bang Nucleosynthesis 
cUb = 0.0214 ±0.0020 |Q| it is clear that we require a dark matter component to explain the 
data. 

Other cosmological probes such as clusters of galaxies and CMB could further constrain 
these models. However such an analysis is beyond the scope of this Letter since it requires 
a detailed re- calculation of, e.g. cluster potentials and CMB perturbations for the models 
discussed here. 

Summary: We have studied the viability of a geometrical explanation for the present 
acceleration of the Universe. This is possible if the Einstein-Hilbert action is supplemented 
with new terms that are negligibly small at high cosmological curvatures but become relevant 
when the curvature of the Universe gets smaller. Despite the phenomenological problems of 
the simplest models, it has been shown that there exists a broad class of modifications of 



8 



gravity that are phenomenologically viable and have accelerated attractors at late-times. In 
this Letter we have performed a detailed numerical analysis of the dynamics of these models. 
We emphasize that this hard numerical problem has not been solved previously. The result 
of this analysis allowed us to compare inverse curvature gravity with Supernovae data. We 
found that SNe data can be fitted in our model without the need of any dark energy and 
getting meaningful constraints in the free parameters. We further have shown that these 
models still require a dark matter component. Of course this latter conclusion does not need 
to hold for more general models, for instance those with n ^ 1. We are planning to study 
more general models and their implications for dark matter in the near future. However, 
we would like to emphasize the generality of our study. We have parameterized all models 
governed by Eq. ([Q) with n = 1. Finally we are currently extending this analysis to CMB 
and cluster datasets, a non-trivial task. This will further constrain these models and maybe 
even distinguish them from dark energy. 
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